Ultracold bosons in a synthetic periodic magnetic field: Mott phases and re-entrant 

superfiuid-insulator transitions 
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We study Mott phases and superfluid-insulator (SI) transitions of ultracold bosonic atoms in a 
two-dimensional square optical lattice at commensurate filling and in the presence of a synthetic 
periodic vector potential characterized by a strength p and a period I = qa, where q is an integer 
and a is the lattice spacing. We show that the Schrodinger equation for the non-interacting bosons 
in the presence of such a periodic vector potential can be reduced to an one-dimensional Harper- 
like equation which yields q energy bands. The lowest of these bands have either single or double 
minima whose position within the magnetic Brillouin zone can be tuned by varying p for a given q. 
Using these energies and a strong-coupling expansion technique, we compute the phase diagram of 
these bosons in the presence of a deep optical lattice. We chart out the p and q dependence of the 
momentum distribution of the bosons in the Mott phases near the SI transitions and demonstrate 
that the bosons exhibit several re-entrant field-induced SI transitions for any fixed period q. We 
also predict that the superfluid density of the resultant superfluid state near such a SI transition 
has a periodicity q (g/2) in real space for odd (even) q and suggest experiments to test our theory. 

PACS numbers: 03.75.Lm, 05.30.Jp, 05.30.Rt 



I. INTRODUCTION 

Several experiments on ultracold trapped atomic gases 
have opened a new window onto the phases of quan- 
tum matter—. A gas of bosonic atoms in an optical or 
magnetic trap has been reversibly tuned between super- 
fluid and insulating ground states by varying the strength 
of a periodic potential produced by standing optical 
wavea^^. This transition has been explained on the basis 
of the Bose-Hubbard model with on-site repulsive inter- 
actions and hopping between nearest neighboring sites of 
the latticci^"— . In fact, experiments on the superfluid- 
insulator (SI) transitions of such bosonic atoms in two- 
dimensional (2D) optical lattices^ is found to agree with 
predictions of theoretical studies of the Bose-Hubbard 
model quite accuratelj«^i^. 

More recently, several experiments have successfully 
generated time- or space- dependent effective vector po- 
tentials for neutral bosons. Such synthetic vector po- 
tentials are created by generating temporally or spatially 
dependent optical coupling between the internal states of 
these bosonic atomsiSr— . We note that this experimen- 
tal technique involves production of a specific effective 
vector potential for the atoms and hence corresponds to 
a fixed gauge. In the simplest experimental setup, these 
vector potentials are typically chosen to represent a con- 
stant magnetic field in the asymmetric gauge. However, 
a few experiments have also generated vector potentials 
which correspond to spatially varying synthetic magnetic 
fieldsi^i. Several theoretical studies have been carried on 
the properties of the bosons in deep optical lattice in the 
presence of a constant synthetic magnetic field^S,. In par- 
ticular, the SI phase boundary has been computed both 
using mean-field theory^i^ and excitation energy calcula- 
tion which relies on a perturbative expansion in the hop- 
ping parameter—. More recently, experimentally relevant 



issues, such as the momentum distribution of the bosons 
in the Mott phase, the critical theory of the SI transition, 
and the nature of the superfluid ground states and collec- 
tive modes near criticality have also been addresse d^^'^^ . 
However, in spite of the possibility of direct experimental 
realizationi^, the phase diagram of these bosons in the 
presence of a spatially dependent magnetic field has not 
been theoretically investigated. 

In this work, we present a theory of the SI transition 
for ultracold bosons in a 2D square optical lattice with 
commensurate filling tiq and in the presence of a periodic 
synthetic vector potential given by A* = (0,^*) with 
A* — AJ sin(27ra;/Z), where I = qa is the period of the 
vector potential, q is an integer, a is the lattice spacing, 
and Aq is the maximum value of the vector potential 
on any lattice site. At the outset, we introduce a di- 
mensionless number p = 2TTq*AQa/hc, (where q* is the 
effective charge of the bosonsii, c is the speed of light, 
and h = 2Trh is the Planck's constant) which will be 
used in the rest of this work to characterize the strength 
of the vector potential. We first consider the problem 
of non-interacting bosons in a lattice in the presence of 
such a periodic vector potential and show that the cor- 
responding single particle Schrodinger equation can be 
reduced to a one-dimensional Harper-like equatio n^^'^^ . 
The solution of this equation yields an energy spectra 
with q bands (with energies e^(k;p) for a = 0..q — 1) 
all of which have a periodicity of 27r/g along k^. The 
lowest of these bands eo(k;p) has, depending onp, either 
a single minimum at k = {kx,ky) = (0,0) or (0,7r) or 
doubly degenerate minima either at (0,0) and (0,7r) or 
at (0, zbfc,™'") where A:™'" can vary continuously as a func- 
tion of p for a given q. The minimum energy of the lowest 
band, emin, turns out to be a non-monotonic function of 
p for a fixed q. Using these properties of the single parti- 
cle energy bands and a strong coupling expansion^ii^, we 



analyze the Mott phase and SI phase transition of these 
bosons in the presence of a deep optical lattice. We show 
that, depending on p and g, the momentum distribution 
of these bosons in the Mott phase near the SI transition 
will exhibit single (double) precursor peak(s) at the posi- 
tion of the minimum (minima) of eQ(k;p). We determine 
the SI phase boundary and demonstrate that the bosons 
exhibit a series of re-entrant field-induced SI transitions 
as a function of the vector potential strength p for any pe- 
riod q. We also construct an effective Landau- Ginzburg 
action for the SI transition and show, by analyzing this 
action at a mean-field level, that the resultant superfluid 
state has a q {q/2) periodic structure in real space for 
any odd (even) q. We show that the reason for such a 
period-halving of the superfluid density for even q can be 
traced back to the properties of the Harper-like equation 
obeyed by the non-interacting bosons. We discuss several 
experiments that can probe our theory. 

The rest of the paper is organized as follows. In Sec. 
ITTl we introduce the relevant tight-binding Hamiltonian 
of the bosons in an optical lattice in the presence of the 
periodic vector potential and obtain the energy spectrum 
when the interaction between these bosons is set to zero. 
This is followed by Sec. lIIH where we introduce the strong 
coupling expansion for the bosons and use it to compute 
the boson momentum distribution in the Mott phase and 
the SI phase boundary. In Sec. IIV) we show that the 
superfluid state into which the transition takes place ex- 
hibits a g-pcriodic superfiuid density. We conclude with 
a discussion of possible experiments to test our theory in 
Sec.lVl 



interacting boson Hamiltonian becomes 

m,n 

+h.c] , 



(3) 



where we have used p = 2Trq*AQa/hc. To obtain the 
spectrum for Hq, we use the identity 



= J2 J.(^)e"^ 



(4) 
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where Jr{z) denotes Bessel functions with integer r, and 
write Jio in momentum-space representation as 

V-a = -t^y^ \2cos{k^)b^{k^,ky)b{kx,ky) 

k 

oo 

+ ( Y. Mp)e-"'''b\k,,ky)b{k^ + 2TTr/q,ky) 
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+h.c 

= -t 



' Y, [2(cos(fc:,) + So{p) cos{ky)) 
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xb\kx,ky)b{kx,ky) + y^5'r(p)e' 

r=l 

xb\k.^, ky)b{k^ + 2nr/q, ky) + h.^ 






(5) 



where b{kx,ky) = ^j^. exp(i(fca;m -|- kyn))bmn- In Eq. [SJ 
Sr{p) is given by 



II. NON-INTERACTING BOSON SPECTRUM 



Sr{p) = Y '^l^+riP)' 



(6) 



The Hamiltonian of a system of bosons in the presence 
of an optical lattice and a synthetic periodic vector field 
is given by^'^'^'i^'i^ 

n = Ytrr'blbr'+Y^-Hnr + -Mnr-l)] (1) 



where /i is the chemical potential, U is the on-site Hub- 
bard interaction, bj. {fir — 6j6r) is the boson annihilation 
(density) operator, the hopping matrix t'^.^, is given by 



= -t'e~"i'^'^ 



A* -dl/hc 



(2) 



if r = (a;,y) = {m,n)a and r' are nearest neighboring 
sites and is zero otherwise, and t' is the hopping ampli- 
tude of the bosons between the nearest neighboring sites. 
In the rest of this work, we set the lattice spacing a, S, 
and c to unity. Our aim in this work is to analyze the 
phases of %. 

To this end, we first analyze the boson spectrum in 
the non-interacting limit /i = C7 = 0. In this case, non- 



where n takes integer values, Srip) = Sr+q{p), and we 
have used the 27r periodicity oib{kx, ky): b{kx:-\-2'K, ky) — 
b{kx,ky). Note that for even q, Sr = for all odd inte- 
ger r which follows from the well-known property of the 
Bessel functions Jn{p) = {—^)^J-n{p) for any integer n. 
The Schrodinger equation obtained from Eq. [1] can be 
written by expressing the eigenfunctions asi^ 
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= ^Va6^(fc. + 27ra/g,fc.y)|0), 



(7) 



Q = 



where ipa = "fpa+q, and obtaining the equations of tpa 
from T-Lolip) — E\^p). This yields a one-dimensional 
Harper-like equation for ipa 

e^Ja = -t'\2{cos{tj:) + So{p)coa{ky))ipa 

q-1 

+ Yl ^r{p) (e-^'^" Va+r + e*'^" Vc-r) ] ■ (8) 
r=l 

Eq. [5] can easily be cast in the form of g x g dimen- 
sional Hermitian matrix equation A'(k;p)'(/; — eip. The 



eo''(k;l) 




FIG. 1: (Color online) Plot of ej, (k; 1) showing a single min- 
ima at {kx, ky) = (0, 0). 



diagonal elements of A^(k;p) are given by A^„(k;p) = 
—2{cos{kx + 27r(ri — 1)/?) + S'o(p) cos(fcy)) and the off- 
diagonal elements by A^,_„+^(k;p) = A^+^„(k;p) = 
—Sr{p)e~^^'« . The difference of A''(k;p) with its counter- 
part in the constant magnetic fieldi^ is two-fold. First, 
A'?(k;p) no longer remains a tri-diagonal matrix. How- 
ever, the 27: / q periodicity of its eigenvalues, which is a 
consequence of the periodicity of the magnetic field, is 
still retained. This property is most easily seen by not- 
ing that a shift of fc^ — > fcx + 27r/q in Eq. [7] amounts 
to a shift of ^Ja -^ V'a+i- Second, for even g, where 
A^^„+,(k;p) = Ar+,,„(k;p) = for aU odd r, Ank;p) 
separates into two block-diagonal matrices of dimension 
q/2 leading to q/2 non-zero elements of the eigenvector 
-0 for any eigenvalue e. Note that for q = 2, which cor- 
respond to A* = on all sites, we have So{p) = 1 and 
Si (p) = so that Eq. [5] reduces to the standard tight- 
binding Hamiltonian in zero magnetic field. 

For q > 3, a straightforward numerical diagonaliza- 
tion of A leads to q energy bands with energy dispersions 
e^(k;p), where a = 0..q — l, which have a period of 2Tr/q 
along kx ■ This periodicity is a manifestation of the q-fold 
folding of the Brillouin zone due to the presence of the pe- 
riodic vector potential. The lowest energy band eg (k; p) , 
shown in Fig. [1] for p = I and q ~ 3, displays a single 
minima at (kx, ky) — (0,0) within the magnetic Brillouin 
zone {—■n/q <kx< Tr/g and — tt < ky < tt). This minima 
structure changes with increasing p as shown in Fig.[2]for 
g = 3, 4, 5 and 6. For g = 3, 5 and 6, we find that beyond 
a critical strength of the vector potential pi{q), eo(k;p) 
has two minima at the (0, ±fc™'"(p)). As p is increased, 



fc™'" increases monotonically from to tt until it reaches 
t: at p = p2 (q) , where we recover the single minima struc- 
ture of eQ(k;p) with the minima at (0, tt). As p is further 
increased, till a value P3{q), fc™'" remains at n. Beyond 




FIG. 2: (Color online) Plot of fc™'° (red solid line) and 
Emin/i' (black dashed line) as a function of the vector poten- 
tial strength p for q — 3 (left top panel), 4 (right top panel), 
5 (left bottom panel), and 6 (right bottom panel), showing 
non-monotonic periodic behavior. See text for details. 



Psiq), for g = 3 and 6, we find that eQ(k,p) again has 
two minima at (0,±A;™™(p)) and fc™'"(p) monotonically 
decreases from tt to as p is increased. For g = 5, be- 
yond ps^q), we find a discontinuous change in A:™'" from 
TT to 0, and e^Ck^p) retains its single minima structure. 
For q = 4, we always have a single minima of eQ(k;p) 
at {kx,ky) = (0,0), except at p = nn where there are 
two degenerate minima at (0, 0) and (0, tt). We also note 
from Fig. [21 that the minimum value of the energy, emin, 
is a non-monotonic function of p for all q < 6. We have 
checked that these features remain qualitatively similar 
for q > 6 and we shall not discuss those cases further 
here. In the next section, we shall utilize these proper- 
ties of Eq (k; p) to understand the phase diagram of these 
bosons in the presence of a deep optical lattice. 

Before ending this section, we note that there is an 
alternative method of finding the energy eigenvalues of 
the Hamiltonian Eq. [3] by constructing the Schrodinger 
equation in real space and using the q periodicity of the 
eigenfunctions along x. This has been carried out in Ref. 
|20| and yields identical results to the method elaborated 
here. We also point out that, although we have, keep- 
ing in mind the simplicity of experimental realization, 
considered a relatively simple sinusoidal form of the vec- 
tor potential, our method can be easily generalized to 
treat more complicated periodic vector potentials. Also, 
we note that since the vector potential A* is not a gauge 
field, there is no gauge freedom in the choice of the eigen- 
functions (Eq. [7]). Thus the flux of the vector potential 
appears only in the coefficients of the Harper equation 
and not in the choice of the eigenfunctions which is in 
contrast to the case of periodic magnetic fields with gauge 
freedom treated in Ref. uU 




n(k) 




FIG. 3: (Color onlme)Plot of n{kx — 0, ky) as a function of 
ky for different representative value of p, and for q — 3 (panel 
1) and Q = 5 (panel 2) at n/U = 0.414 and t'{p)/t'^{p) = 0.95. 



FIG. 4: (Color online)Plot of n(k) for q = 3 and p = 1 at 
fi/U = 0.414 and t' /t'^ = .95 showing peaks at minima (0, 0). 



III. STRONG COUPLING EXPANSION 

In this section, we analyze the phases of 7i in the Umit 
of t' /U <C 1, where the bosons are in a Mott insulating 
state. We note that the effect of the magnetic field man- 
ifests itself in the first term of Eq. [T] and thus vanishes 
in the local limit {t' = 0). In this limit the boson Green 
function can be exactly computed^i^iii^ and is given, at 
T = 0, by 



Go(iuj„ 



(no + 1) 

itUn — Er, 



no 



Eh 



(9) 



Here w„ denote bosonic Matsubara frequencies and Eh — 
/i — U{no — 1)(-Ep = — M + Uno) are the energy cost of 
adding a hole (particle) to the Mott state. To address 
the effects of the hopping term, we resort to the coher- 
ent state path integral description of these bosons. The 
partition function of the system can then be written as 



Sa = drJ2[rrir)drM^ 



fJ-nrir) 



+ -^nr(T-)(nr(r) -1) 

^1 = / dr^C'CW^r'(T). 

Jo I 



(10) 



Here r is the imaginary time, ip denote boson fields in 
the path integral representation, nr{T) = 'iP*(t)iPj.(t), 
P = l/ksT is the inverse temperature (T), and fc^ is the 
Boltzman constant. Following Ref. 0, we then decouple 
the hopping term introducing a Hubbard-Stratonovitch 



field (/)r (■'■)■ The partition function can then be written 
as 



Z = DiPDi;*D(j)D(j)*e-'-^''+^'^+^'^\ 

S'l = f dTj2{rrir)Mr)+h.c), 
Jo 

S'2 ^ - I rfr^C'V:(r)</>r'(r). 



(11) 



Finally, we introduce a second Hubbard-Stratonovitch 
field i!r{T) and decouple S'2 to obtain 



Z 



Di}Di}*D(j)D(j)*D%l;D%i;*e-'^^°+^''+^'^\ 



S'i = I dTV4,^;(r)^r'(T). 
Jo I 



(12) 



Note that integrating out 4>r [t) in Eq. [T2] would lead to 
the constraint "^r — V'r on Z . It can also be shown that ip 
and f/' fields have identical correlation functions^. Next, 
we follow Refs. .7.. 16. to integrate out the '0 and </> fields 
and obtain an effective action in terms of ^j. The details 
of this procedure has been elaborated in Ref. 0. The 
effective action so obtained is given by^i^^ 

5cff — Sq + Si 
So = [ rq{^^n,'k)[-Go\^LOn)I + A'{^k)]^l;g{iu„,^k), 

Si = g/2 I dr f d\\rg{r,r)^g{r,T)\\ (13) 



n(k,=0,k„) 6 




n(kx=0,k„) 




FIG. 5: (Color online) Plot of n{kx — 0,ky) as a function of 
3 at fi/U = 0.414 and t' /U = 0.04. 



FIG. 6: (Color online) Same as Fig. [5] for q — A. 



where ipq 



= {lpo{kx,ky)..1pq-l{ka:,ky))'^ wlth 

ipa(kx,ky) = ip{kx + 27ra/g, fcy) denoting the q- 
component of the auxiliary field "0 in momentum space, 
/^ = [1/ P) Yl,ur, I d^k/{2'KY, I denotes the unit matrix, 
and g > is the static limit of the exact two-particle 
vertex function of the bosons in the local limit which 
has been computed in Ref. 0. Note that 5*0 reproduces 
exact bosons propagator both in the local (i' = 0) 
and the non-interacting {U = 0) limits and therefore 
provides a suitable starting point for the strong coupling 
approximation. In the next subsection, we shall compute 
the momentum distribution function of the bosons from 

^0. 



A. Momentum distribution of the bosons in the 
Mott phase 

The momentum distribution of the bosons in the Mott 
phase can be computed from Sij^-^ 

n(k) = -lim(l//?)VTrG(zc^„,k), 



G{iuJnM) 



-Go^(z^„)^ + Ank;p)]-i. (14) 



To compute n(k), we note that Gq is independent 
of momenta. Hence finding G(iw„,k) amounts to in- 
verting A'(k;p). To this end we introduce an unitary 
transformation where the transformation matrix [/^(k) 
diagonalizes A''(k;p) to obtain a diagonal Green func- 
tion G'^{iuJn, k) = C/~^(k)G(zw„, k)C/g(k) whose diagonal 
elements are given by 

iu!n + fi + U 



GJ^„(iw„,k) 



-Gq ^{iuj„ 



{iiOn-E^+ik;p)){ioJr,-E^-ik;p))'' 



where we have used the expression of Gq from Eq. [9] and 
E^^ (k; p) denote the location of the poles of the inter- 
acting boson Green function and are given by 

E^^ik;p) = -A* + C/(no-l/2) + e^(k;p)/2±i 

e« (k;p)2 + 4e^(k;p)C/(no + 1/2) + C/2. 

(16) 

Note that i?"^(k;p) can be directly computed from 
the knowledge of the non-interacting boson spectrum 
e^(k;p) derived in Sec. |lTl In particular, the minima 
-E" (k;p) occur in the same position in the magnetic 
Brillouin zone as e^(k;p). Also, as noted in Ref. 0, the 
Mott gap E"~^(k;p) — E"^(k;p) vanishes at the position 
of the minima of e^ (k; p) in the magnetic Brillouin zone 
provided we are at the tip of the Mott lobe where the SI 
transition takes place at constant density. 

The momentum distribution can now be computed as 
n(k) = -limT^o(l//3)Ea; TrG''(ia;„,k) and is given 
bj*l£ 



g-l 



n(k) - J2 



E^-{k;p)+fi + U 






(k;p) 



E^- 



(k;p) 



(17) 



Eg. [T71 shows that the peaks of n(k) occur when the Mott 
gap E"~^{]s.) — E"~{'k) becomes small near the minima 
of e^ (k; p) as the SI transition is approached through the 
tip of the Mott lobe. The minima structure of the non- 
interacting bosons is therefore expected to be reflected in 
the peaks of the momentum distribution of the bosons in 
the Mott phase. In Fig. [H we show a representative plot 
of n(k) as a function of k for q = 3 and p = 1. We find 
that the central peak of the momentum distribution lies 
at (0, 0) in accordance with the position of the minima 

(3) 

(15) of Eq (k, 1). Next, keeping in mind that the position of 
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n(kx = 0,ky) g' 





FIG. 7: (Color onlme)Sanie as Fig. [5]for q — 5. 



FIG. 8: (Color online)The MI-SF phase boundary for g = 3 
and p — 1. 



the minima of e^(k;p) always occur at fc^; = 0, we plot 
the momentum distribution n{kx — 0,ky) as a function 



of ky (for fixed t'{p)/t'^{p) = 0.95 and 9 = 3,5) for several 
representative values of p in Fig. |3l Fig. |3] clearly shows 
that as p increases, the peak structure of the momentum 
distribution changes from a single peak at fcj, = to two 
split peaks at ky = ±A:™'"(p) and finally to a single peak 
a.t ky = TT. Finally in Figs. [SI [5] and [71 we plot n{kx = 
0, ky) for q — 3, 4, and 5, as a function of ky and p for a 
fixed t' = 0.04C/. Note that for these plots, the proximity 
of the system to the tip of the Mott lobe changes with 
p since t'^ is a function of p. These plots again reveal 
the change in the peak structure of n{kx: ~ 0, ky) as a 
function of p. 



B. Re-entrant SI transitions 

The critical hopping t'^ for the MI-SF transition as a 
function of /i can be determined from the condition^^ 



rq{p) 



-Go-i(z^„ = 0) + 6f„i„(p) = 0. (18) 



The SI phase boundary so obtained is shown in Fig. [T] 
for g = 3 and p = 1 in Fig. [8| and displays the usual 
Mott lobes. The difference of the present case here with 
the SI transitions studied earlier— iiiSiii^ arises due to the 
non-monotonic p dependence of e^in(p). This point is 
demonstrated in Fig. [9| for q = 3, 4, 5, and 6 by plotting 
t'^{p) as a function of p for hq — 1 and ^ = /xtip. We find 
that t'^{p) is a non-monotonic function of p and t'^{p) > 
t'^{0) for all p, Consequently, varying p at a fixed value of 
t' > t'^{0) leads to a series of field-induced re-entrant SI 
transitions for any q. This is schematically marked by the 
red-dotted line in Fig. [9] We note that such re-entrant 
transitions as a function of the magnetic field strength 
are not present for SI transitions in a constant magnetic 
fieldl^. 



IV. THE SUPERFLUID PHASE 

At t' — t'^{p), it becomes energetically favorable to cre- 
ate particles/holes at the minima of the energy dispersion 
of the bosons leading to the destabilization of the Mott 
phase. The Landau- Ginzburg theory of the resultant su- 
perfluid phase can be expressed by long-wavelength bo- 
son fields around these minima. In the present case, there 
are either one or two degenerate minima of the boson en- 
ergy spectrum in the magnetic Brillouin zone leading to 
a Landau-Ginzburg theory of one or two low-energy bo- 
son fields^'i'Sii^. We shall first consider the case with a 
single minima either at (0, 0) or (0, n) which occurs for 
specific ranges of p for all q as discussed in Sec. [Ill In 
either case, the boson field can be written as 

ip{r,t) = xo{r;p)(p{r,t), 

"9-1 



xo(r;p) 



J2^'='^P^^ 



27riax/q 



e''=r"V(r,i), (19) 



where ipaip) denotes the components of eigenvectors of 
Aq(k;p) at k^ — 0, ky = fc™™ which can be either or 
TT for a fixed p, and xo(r;p) denotes the corresponding 
wavefunction in real space. Thus the superfluid density 
can be written as 



ps{r) = m\'^ 



9-1 

^Vc.(p)e""""/« 

Q = 



l^oP, (20) 



where tpo — {(p{r,t)) ^ for t' > t'^{p). Note that ps 
is independent of y irrespective of the value of fc™™, but 
displays spatial variation along x. Further, as discussed 
in Sec. [ill for even g, only q/2 of the components tpa 
(corresponding to either even or odd integers a) will be 
non-zero. Consequently, we expect the period of Ps{x) to 
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FIG. 9: (Color online) Plot of the critical hopping strength 
t'c{p) as a function of p at the tip of the Mott lobe (/i — fitip = 
0.414f/) for (J = 3 (left top panel), 4 (right top panel), 5 (left 
bottom panel), and 6 (right bottom panel). The red-dashed 
line is a guide to the eye showing reentrant SI transitions as 
p is varied at fixed t' > t'c{p = 0). 




FIG. 10: (Color online) Plot of the superfluid density 
Ps{x)/ps{0) as a function of x for g = 3, 5 (left panels) and 
q = 4, 6 (right panels). Note that the superfluid density dis- 
plays a q periodic pattern for odd qs and a q/2 periodic pat- 
tern for even qs. p is set to 0.5 for all plots. 



be halved. A plot of the renormalized superfluid density 
Ps{x)/ps{0), plotted in Fig. [10] for p = 0.5 and g = 3, 4, 
5, and 6, confirms this expectation. The presence of the 
periodic vector potential leads to a q- periodic pattern 
with q — 2 small and one large peak in the superfluid 
density along x for all odd q as shown in the left panels 
of Fig. [ini In contrast, the superfluid density for even 
q displays a q/2 periodic pattern. Note that this period 
halving leads to identical superfluid density patterns for 
vector potentials with periods q and 2g for all odd q. This 
feature is clearly demonstrated in the top left {q — 3) and 
the bottom right {q = 6) panels of Fig. [TUl 

Next, we derive the effective low-energy Landau- 
Ginzburg theory. To this end, we substitute Eq.[Tn]in Eq. 
1131 and obtain the effective low-energy Landau-Ginzburg 
action in terms of the ip fields. The details of this pro- 
cedure is charted out in Ref. [la- The resultant action is 
given by 



S', 



LG 



(frdt 



ip*ir,t)[Kodf + iKidt + rg{p) 



-v,{pndl + d^^)]^ir,t) + ^Mr,t)\^ 



,(21) 



where Ko = l/2d^Ga^/duj\=o = "0(^0 + l)U^/{^^ + 
Uf, K, = dGo'/dio\^=o = 1 - noino + 1)C/V(m + U)\ 
and Vq{p)'^ = V^emin(k;p)/2, r,(p) is given by Eq. [JHl 

and g' = gEl~y=o \xo{r;p)\'^/q^. At the tip of the Mott 

lobe, where p = ^tip = U{^/r^i{T^i'+l) - 1), Ki = 0. 
Thus we have a critical theory with dynamical critical 
exponent z — 1. Away from the tip, Ki ^ rendering 
z = 2. Thus the critical theory turns out to have similar 
exponent as in the case without magnetic field^. 



Finally, we briefly comment on the case where there 
are two degenerate minima either at (0, ±fc""") or at 
(0,0) and (0,7r). In this case, i^{r,t) — Xo {tc; p)'p+{r ,t) + 
Xo^(r;p)v?-(r,i) where Xo C'^) denotes the eigenfunctions 
of A(k; p) in real space at (0, ±fc™'") and ip± (r, t) denotes 
low-energy fluctuating fields about the minima. Substi- 
tuting this expression of ^ in Eq. I13[ and following the 
coarse-graining procedure detailed in Ref. Ila . we find 
that for all q and p, the superfiuid phase corresponds 
to the condensation of only one of the low-energy fields: 
(ip+) = 0, (ip-) ^ OT (^_) = 0, {ip+) ^ 0. Thus the 
effective Landau-Ginzburg action in these cases is quali- 
tatively similar to Eq. [211 The superfluid density, plotted 
in Fig. [11] for q = 3,5 and p = 2.5, shows similar q peri- 
odic pattern as observed in Fig. [10] for odd q. 



V. DISCUSSION 

There are several possible experimental verifications 
of our theory. First, we suggest measurement of n(k) 
for the bosons in the Mott phase near the transition as 
done earlier in Ref. [^ for 2D optical lattices without the 
synthetic magnetic fleld. Our prediction is that the peak 
structure of the momentum distribution along fc^; = at 
a fixed t' /U near t'^ would be similar to those shown in 
Figs. [5] [7] In particular the shift in the peak position 
of n(0,ky) with p and change from a single to double 
peak structure as a function of p should be observable in 
such experiments. Second, the re-entrant SI transition 
can also be verified by measuring n(k) as a function of 
p by fixing t' > t'^{p = 0) as shown in Fig. [5] Finally, 
the spatial variation of the superfluid density can also be 




observed by measuring n(k) in the superfluid phase. 



In conchision, we have analyzed the MI-SF transition 
of uhracold bosons in a 2D optical lattice in the presence 
of a synthetic periodic magnetic field. We have shown 
that the precursor peaks of the momentum distribution 
in the Mott phases can be tuned by the strength p of 
the synthetic field. We have also demonstrated that the 
bosons, in the presence of such a periodic synthetic mag- 
netic field, show a series of field-induced re-entrant SI 
transitions, and that the superfluid density in the SF 
phase near criticality shows q {q/2) periodic spatial pat- 
tern for odd (even) q. We have suggested several experi- 
ments which can test our theory. 



FIG. 11: (Color online) Same plot as Fig. [TOl for o = 3 
{kf'' = 2.2)) (left panel) and <j = 5 {kf" = 1.89) (right 
panel) showing q periodic patterns, p = 2.5 for both cases. 



K.S. thanks R. Shankar for discussions and DST, India 
for flnancial support under Project No. SR/S2/CMP- 
001/2009. 



^ M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002); 
^ C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, and 

M. A. Kasevich, Science 291, 2386 (2001). 
^ M. P. A. Fisher, P. W. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 

* S. Sachdev, Quantum Phase transitions, Cambridge Uni- 
versity Press, (1999). 

^ D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. 

ZoUer , Phys. Rev. Lett. 81, 3108 (1998). 
^ K. Seshadri, H. R. Krishnamurthy, R. Pandit, and T. V. 

Ramakrishnan, Europhys. Lett. 22, 257 (1993); 
■^ K. Sengupta and N. Dupuis, Phys. Rev. A 71, 033629 

(2005). 

* I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. Rev. 
Lett. 98, 080404 (2007). 

® J. Freericks, H. R. Krishnamurthy, Y. Kato, N. 
Kawashima, and N. Trivedi, Phys. Rev. A 79, 053631 
(2009). 

^° D. Jaksch and P. ZoUer, New J. Phys. 5, 56 (2003); E. 
Mueller, Phys. Rev. A 70, 041603(R) (2004); K. Osterloh, 
M. Baig, L. Santos, P. ZoUer, and M. Lewenstein, Phys. 
Rev. Lett. 95, 010403 (2005); N. Goldman, A. Kubasiak, 
P. Gaspard, and M. Lewenstein, Phys. Rev. A 79, 023624 
(2009); I. B. Spielman, Phys. Rev. A 79, 063613 (2009). 

" Y-J. Lin, R. L. Compton, A. R. Perry, W.D. Phillips, J.V. 



Porto, and I. B. Spielman, Phys. Rev. Lett. 102, 130401 

(2009). 

Y.-J. Lin,R. L. Compton, K. Jimenez-Garca, J. V. Porto, 

and I. B. Spielman, Nature 462, 628-632 (2009). 

N. Goldman, A. Kubasiak, A. Bermudez, P. Gaspard, M. 

Lewenstein, and M. A. Martin-Delgado, Phys. Rev. Lett. 

103, 035301 (2009); I. Satija, D. C. Dakin and C. W. Clark, 

Phys Rev lett, 97, 216401, (2006); S-L Zhu, Hao Fu, C- 

J. Wu, S.-C. Zhang, and L.-M. Duan, Phys. Rev. Lett. 

97, 240401 (2006); H. Zhai, R. O. Umucalilar, and M. O. 

Oktel, Phys. Rev. Lett. , 104, 145301 (2010). 

R. O. Umucahlar and M. O. Oktel, Phys. Rev. A 76, 

055601 (2007); E. Lundh, EuroPhys. Lett. 84, 10007 

(2008). 

M. Niemeyer, J. K. Freericks, H. Monien, Phys. Rev. B 60, 

2357 (1999). 

S. Sinha and K. Sengupta, arXiv: 1003:0258 (unpublished). 

S. Powel, R. Barnett, R. Sensarma, S. D. Sarma, 

arXiv: 1004:0701 (unpublished). 

D. Hofstadter, Phys. Rev. B 14, 2239 (1976). 

M. Kohmoto, Phys. Rev. B 39, 11943 (1989). 

G-Y Oh, Phys. Rev. B 60, 1939, (1999). 

F. Guil, Journal of Math. Phys. 48, 033503 (2007); B.A. 

Duvrovin and S.P. Novikov, Sov. Phys. JETP 52, 511 

(1980). 



